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We extend our maximum likelihood method for reconstructing the cluster-mass cross-correlation from cosmic 
microwave background (CMB) temperature anisotropies and develop new estimators that utilize six different 
quadratic combinations of CMB temperature and polarization fields. Our maximum likelihood estimators are 
constructed with delensed CMB temperature and polarization fields by using an assumed model of the conver- 
gence field and they can be iteratively applied to a set of clusters, approaching to the optimal condition for the 
lensing reconstruction as the assumed initial model is refined. Using smoothed particle hydrodynamics simu- 
lations, we create a catalog of realistic clusters obtainable from the current Sunyaev-Zel'dovich (SZ) surveys, 
and we demonstrate the ability of the maximum likelihood estimators to reconstruct the cluster-mass cross- 
correlation from the massive clusters. The iTT temperature estimator provides a signal-to-noise ratio of a factor 
3 larger than the iEB polarization estimator, unless the detector noise for measuring polarization anisotropies is 
controlled under 3/^K. 
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I. INTRODUCTION 

In the past few years we have seen rapid development in 
the measurements of cosmic microwave background (CMB) 
anisotropies on arcminute scales. Higher precision measure- 
ments of CMB temperature anisotropies are available from 
the completed missions such as the Arcminute Cosmology 
Bolometer Array Receiver (111 and the Cosmic Background 
Imager [2] and the ongoing experiments with better angu- 
lar resolution and multi-frequency channels such as the At- 
acama Cosmology Telescope (ACT, yj]) and the South Pole 
Telescope (SPT, |4]). While most of the current polariza- 
tion experiments are aimed at measuring CMB polarization 
anisotropies on large scales, ACT and SPT will be capable 
of measuring CMB polarization anisotropies on arcminute 
scales. 

Especially, ACT and SPT are designed for measuring the 
Sunyaev-Zel'dovich (SZ) effect arising from the Thomson 
scattering off hot electrons in clusters and thereby detecting 
massive clusters. These SZ surveys can explore the growth of 
structure by measuring the abundance of massive clusters and 
the expansion history of the universe by measuring their cor- 
relation functions to probe the evolution of dark energy (see, 
e-g-^ |5il). With precisely measured distance and their om- 
nipresence behind all the clusters, CMB anisotropies can be 
used as a distant background source for weak lensing mea- 
surements, providing the cluster-mass cross-correlation and 
complementing the measurements of cluster mass and their 
abundance from the same SZ surveys. Moreover, weak lens- 
ing measurements of the CMB can be combined with galaxy 
weak lensing measurements of the same clusters from opti- 
cal follow-up surveys to measure the source distance ratios in 
a model-independent way, constraining the expansion history 
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of the universe ia] . 

Weak lensing of the CMB by large-scale structure is effi- 
ciently probed by using quadratic estimators, measuring the 
deviation of the correlation function of CMB anisotropies 
from the otherwise statistically isotropic correlation Q . How- 
ever, it was shown that the standard method of using quadratic 
estimators is compromised 18] on cluster scales and needs 
additional free parameters 131 for calibration, as this method 
rests on the linear approximation in the lensing effect which 
breaks down near massive clusters. This problem was tackled 
by delensing the observed CMB temperature field and analyz- 
ing the likelihood of the delensed fields [ 10|] . 

We showed that if the assumed initial model for delensing 
is a good approximation to the true underlying matter distri- 
bution, our new estimator based on the delensed CMB tem- 
perature fields becomes an optimal estimator and it can be 
iteratively applied to given measurements of CMB tempera- 
ture anisotropies, until the assumed initial model converges 
to the true matter distribution and the likelihood is maxi- 
mized. Here we extend our maximum likelihood estimation 
method for reconstructing the cluster-mass cross-correlation 
with CMB temperature anisotropies and apply it to CMB po- 
larization anisotropies, yielding new maximum likelihood es- 
timators with six different quadratic combinations of CMB 
temperature and polarization fields. We demonstrate their ap- 
plicability by using realistic clusters that can be found in the 
current CMB experiments. 

The rest of the paper is organized as follows. We briefly 
review the formalism for weak lensing of CMB anisotropies 
and derive our maximum likelihood estimators in Sec. [Ill In 
Sec. [Ill] we describe our numerical simulations to model re- 
alistic clusters and we construct a catalog of massive clusters 
with CMB secondary anisotropies. We reconstruct the cluster- 
mass cross-correlation by using our maximum likelihood es- 
timators and quantify their signal-to-noise ratios in Sec. lIV Al 
and we investigate the impact of the kSZ contamination in 
Sec. lIV Bl We summarize our findings and conclude in Sec.lV] 



For illustrative purposes, we adopt a flat ACDM universe with 
cosmological parameters (flmh^ = 0.127, fl^h^ — 0.0222, 
h — 0.73, Ug = 0.95, (7h = 0.78), consistent with the recent 
estimation (e.g., llTll|l2|l ). 



11. FORMALISM 

A. Weak lensing of CMB polarization 

Gravitational lensing is a surface brightness conserving 
process and it simply redistributes the intrinsic CMB temper- 
ature and polarization fields. Here we use T, Q, and U to 
represent the CMB temperature and the Stokes parameters for 
the CMB polarization in units of /iK and use tildes to indi- 
cate that the corresponding field is gravitationally lensed. The 
lensed CMB temperature and polarization fields at the angular 
position n on the sky are therefore 

f(n) = T[n + V0(n)j , (1) 

(Q±zC/)(n) = (Q ± zt/) [n + V^(n)l , 



and the projected potential 







dD^^^P{Dh,D) 
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describes the deflection angle V0(fi) on the sky, where the 
gravitational potential is V', the comoving angular diameter 
distance to the CMB last scattering surface is D,t = 14.12 Gpc 
II12I1 . and V is the derivative with respect to unit angular vector 
fi. The projected potential cj) is related to the convergence k 
as V^0 = —2k and it is further related to the physical matter 
density Spm projected along the line-of-sight as 



K(n) 



dD ,, *,„„^ AnG SprniDh, D) 







{i + zyD, 



(3) 



As we are interested in reconstructing the cluster-mass 
cross-correlation, we adopt the flat sky approximation and ex- 
press quantities of interest in Fourier space. To linear order 
in (p, the lensed CMB temperature field can be written as 



Ti = Ti- 



(^ (I2 • li) 0b Ti, , 



(4) 



with I2 — 1 — li and this equation is also valid for the 
lensed Stokes parameters Q\ and U\ in Fourier space. How- 
ever, since density fluctuations excite only curl-free polariza- 
tion, the CMB polarization fields are better described by two 
parity eigenstates, E- and _B-modes lll3|] as (Qi ± iU\) = 
{E\ ± iBi)e^'^^'f^ with the phase {p\ of the wavevector 1. The 
lensed CMB polarization fields are therefore 
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and the lensed CMB power spectra are related to the intrinsic 
CMB power spectra as 



Cf = [I- fR] Cj 



(2^)2 

}-2 
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(6) 



Cf ^ [1 -f R] Cf + J ^{h -hfctC^^ cos' 2A^,,, 
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cF = 
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where A(pij — Lp\^ — ip\ and R — (l/47r) fdlnZ I'^Cf is the 
half of the rms deflection angle (see, e.g., lll4l[l5ll ). 

Finally, we assume the Gaussian random noise for the de- 
tector and the Gaussian beam of the telescope. Therefore, the 
power spectrum of the detector noise is lUol 



C, 






pix ''"pix J 



(7) 



and the observed CMB temperature field and its power spec- 
trum are 






(8) 



where cTpix is the rms noise of the detector, ilpix is the solid 
angle subtending each pixel of the detector, N^ is the Fourier 
mode of the detector noise, and the full-width half-maximum 
(FWHM) of the telescope beam is 6'fvvhm = ctj, v81n2 . 
The observed CMB polarization fields are also described by 
Eqs. ([HJ, while cTpix and cr;, in CMB polarization experiments 
may differ from those in CMB experiments for measuring 
temperature anisotropics. 



B. Improved quadratic estimators 

Matter fluctuations along the line-of-sight deflect CMB 
photons and this process imprints a deviation of the CMB two- 
point statistics from statistical isotropy. Quadratic estimators 
are often used to measure the deviation, but they become pro- 
gressively biased as the lensing effect increases. This problem 
can be overcome by using improved quadratic estimators that 
are constructed by using delensed CMB fields [101 . Here we 
briefly review the standard quadratic estimators 091] using the 
CMB polarization fields, on which our improved quadratic es- 
timators are based, and then we describe how our improved 
quadratic estimators can be used to reconstruct the cluster- 
mass cross-correlation. 

With the CMB temperature and polarization fields, we 
can construct six convergence estimators K^^(n) that use 
two CMB temperature and polarization fields X,Y = 
T, E, B, permitting repetition (hence they are quadratic es- 
timators). The estimators K^^(n) are symmetric with two 
fields X and Y interchanged and we use quantities with a 



hat to indicate that they are estimators thereof (not to be con- 
fused with unit angular vectors on the sky). The six conver- 
gence estimators are uniquely determined with the require- 
ment that each estimator be unbiased (^'^^(n)) = K(n) over 
an ensemble average of the CMB temperature and polarization 
fields X and Y, and the variance of the estimator be minimal 
10], 
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) = (27r)2,5^(l-l')(cr+iVr' ) 



(9) 



These estimators can be easily computed in configuration 
space by using the two Wiener-filtered functions 
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with the phase angle (p\ of the wavevector 1 and two phase 
factors for y = £■, B in the braces. Note that C/^^ == C{^^ 
for Y ~ B and the phase factor in the braces is unity for 
Y ^T. The convergence estimator is then 
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(11) 
The normalization coefficients Ni^^ are related to the 
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noise power spectrum N^" of XY estimators k (fi) as 
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I Nf- /4, and they can be obtained as 
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Wiener-filtered functions Gxy(n) and WV(fi) but with the 
delensed CMB temperature and polarization fields X'^ and Y"^ 
rather than X°^<^ and Y°^^ themselves, and Eq. (fTTT i yields the 
change AK^^(n) in the convergence field with respect to the 
assumed initial model k™ (n). The resulting convergence field 
K™(n) + Ak-^^ (n) then serves as a refined initial model for 
the next iteration and this process is iterated until the numeri- 
cal convergence is achieved AK^^(n)/K'"(n) ^ 1. In prac- 
tice, one can try different initial models for a faster numerical 



convergence, when Ak (n) 



'(n). 



While our improved quadratic estimators take the same 
functional form as the standard quadratic estimators (thereby 
its name is inherited), there exist critical differences: Our 
improved quadratic estimators are based on the maximum 
likelihood of the delensed CMB temperature and polarization 
fields, and the iteration process is indeed the standard Newton- 
Raphson method for maximizing the likelihood of each esti- 
mate given the observed CMB temperature and polarization 
fields. Since the initial model K'"(n) that depends on the 
previous iteration process is refined and the refined model 
K™ + Ak becomes a new initial model in the following it- 
eration, the improved quadratic estimators are in fact ratio- 
nal functions of the CMB temperature and polarization fields, 
rather than quadratic functions, making full use of the infor- 
mation contained in the likelihood. 

Furthermore, our improved quadratic estimators are free 
from the approximation that the lensing effect is weak, the 
breakdown of which plagues the standard quadratic estima- 
tors. Consequently, there is no arbitrary cutoff scale used in 
modified quadratic estimators [9] that breaks the symmetry in 
the convergence estimators when the CMB temperature and 
polarization fields X and Y are interchanged. Hereafter, we 
refer to improved quadratic estimators using delensed CMB 
temperature and polarization fields X and Y as iXY estima- 
tors. 



with 1 = li 
where 

fTT 
Jh.h 



h, Aif = if^, - (^1,, and (Xi^Fi,) = /,fi^ 0i , 



III. MASSIVE CLUSTERS 



{\-h)Cl + {\-h)Cl, 



(13) 



fll = il-h)Cf;cos2A^+il-h)Cf;, 

fll = (l.li)C,^sin2A^, 

/iff. - [(Mi)Qf + (l.l2)Qf]cos2A^, 

/iff. = (Mi)C;fsin2A^ 

(see DSHTtII for details). Due to the vanishing signal-to-noise 
ratio, no quadratic BB estimator is used. 

Improved quadratic estimators are similar in many aspects 
to the above standard quadratic estimators. Our estimation 
process is as follows: Improved quadratic estimators take an 
initial model K™(n) of the convergence field and we first com- 
pute the delensed CMB temperature and polarization fields 
X'^{s) — X°^^{h) by solving the lensing equation s ~ 
n + V(/i™(n), where the initial model i5i'"(fi) for the projected 
potential is related to K™(n) as V20™(n) = —2 K™(n) . Im- 
proved quadratic estimators are constructed by using the same 



We describe our model for massive clusters in Sec. IIIIAI 
and discuss possible contaminants for lensing reconstruction 
arising from secondary anisotropics in the CMB temperature 
and polarization fields in Sec. lIIIBl Tests of the applicability 
of our improved quadratic estimators to reahstic clusters are 
presented in Sec. HV] 



A. Numerical simulation 

Here we model massive clusters using the numerical simu- 
lations of Springel and Hernquist lUSll . We use the largest vol- 
ume simulation (G-series) among a number of simulation runs 
with varying mass and spatial resolution. The smoothed par- 
ticle hydrodynamics (SPH) simulations were performed using 
the parallel GADGET code Iil9il . employing the variational 
formulation of Springel and Hernquist ||20|] . We focus on the 
simulation with 2 x 324^ of gas and dark matter particles 
in a comoving cubic volume (100/i~^Mpc)^. Dark matter 
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FIG. 1 : (color online) Average convergence fields fi:(n) of a 40' x 40' 
region around massive clusters at 2: = 1. Each panel shows K(n) ob- 
tained by averaging 1 (top left), 10 (top right), 100 (bottom left), and 
1000 (bottom right) identical but randomly rotated regions around 
the massive cluster from the smoothed particle hydrodynamics sim- 
ulation, which represents different clusters of the same mass. The 
cluster has of mass M = 3.7 x 10 h"^ Mq and virial radius 
Rvii = 1.9/i~^Mpc. Circles show 7?vir on the sky subtended by 
2. '7 and the dotted line in the upper left panel shows the boundary of 
the cluster members identified by the friends-of-friends algorithm. 
Once many lines-of-sight are stacked, the average convergence field 
restores the spherical symmetry. 



halos are identified by applying the friends-of-friends (FoF) 
group finding algorithm |21] to the dark matter distribution 
with a comoving linking length of 0.2 times the mean inter- 
particle separation 62/i~^kpc. Although the simulation was 
run in a flat ACDM universe with slightly different cosmolog- 
ical parameters {Q,m = 0.3, Vt^ = 0.7, fib = 0.04, Us = 1, 
(Tg = 0.9), we simply retain the physical properties of the 
simulation such as particle mass, position, and velocity. 

To generate a catalog of massive clusters that can be found 
in the current Sunyaev-Zel'dovich surveys, we first select the 
most massive halo in the simulation output at z = 1; The 
massive halo contains ^150,000 dark matter particles, corre- 
sponding to virial mass M = 3.7 x lO^'^h^^AlQ and radius 
i?vii = 1.9/i~^Mpc. Every particle in the simulation is then 
shifted to have the massive halo at the center of the simulation 
box by using the periodic boundary condition, and the whole 
simulation box is randomly rotated to provide the distant ob- 
server with different lines-of-sight to the massive cluster. We 
treat the massive cluster seen at different lines-of-sight as in- 
dependent clusters of the same mass at z = 1. This process 
is necessary to generate a large number of massive clusters 
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FIG. 2: Convergence profiles ti{9) around massive clusters at z — 1. 
Profiles are obtained by averaging the stacked convergence fields 
K(n) in Fig.[T]over the annulus 6 around the cluster center. Thin and 
thick lines show the convergence profiles obtained by counting only 
the cluster members identified by the friends-of-friends algorithm 
and by counting all the particles within the annulus, respectively. 
The thick solid line represents the cluster-mass cross-correlation and 
it deviates substantially from the average cluster mass profile (thin 
solid) around the virial radius. The average cluster mass profile is 
well approximated by a projected NFW profile truncated at the virial 
radius (vertical dotted). Note that dot-dashed lines are largely ob- 
scured by solid lines, indicating fast convergence of n{0). 



using SPH simulations. Since individual clusters identified 
by the FoF algorithm lack spatial symmetry, snapshots of the 
cluster at different lines-of-sight are relatively independent of 
each other as far as the projected matter density is concerned. 
Figure [T] shows the convergence fields K(n) of a 40' x 40' 
region on the sky around the massive cluster with the virial 
radius indicated as the circle at the center. With the fixed 
distance Di, to the background source, the critical (physical) 
surface density is just a function of the redshift of the lens- 
ing cluster l/Scrit = 47rGL'(L'* - D)/Di,{l + z) and the 
convergence is the ratio of the projected mass density S(n) 
to the critical surface density K(n) = I](n)/I]crit ■ The co- 
moving angular diameter distance to the redshift slice z = 1 
IS, D — 2400/i^^Mpc and the corresponding critical surface 
density is Edit = 1800 /iMopc"^. Here we keep fixed the 
detector scale 0.'2 per pixel, smaller than the typical tele- 
scope beam size. Each dark matter particle in the simula- 
tion is of mass mdm = 2.1 x IO^/i^^A/q and it contributes 
K = 6.0 x 10^^ per pixeP at z = 1. Many lines-of-sight 
are randomly generated and each panel shows ^(n) averaged 
over 1,10, 100, and 1000 different lines-of-sight to the mas- 
sive cluster, mimicking the process of stacking many different 
clusters at z — 1. Once many clusters are stacked, the av- 
erage convergence field (bottom right) restores the spherical 
symmetry. 
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FIG. 3: (color online) Kinetic Sunyaev-Zel'dovich (kSZ) effect. 
Top left: Thomson scattering optical depth r of the same region in 
Fig. \T\ We assume the helium mass fraction is Iho ~ 0.24 and 
helium atoms are doubly ionized. Top right: Line-of-sight veloc- 
ity field wios of dark matter particles. Bottom left: The kSZ ef- 
fect ATksz/r = — T fios/c. Electrons are assumed to follow the 
dark matter distribution. Bottom right: ATksz from the kSZ effect 
is superimposed to the intrinsic CMB temperature field. The hori- 
zontal color bar represents the scales of each panel in the range of 
-4.4 < logr < -2.0, -700 kms"^ < wio^ < 560 kms~\ 
-16 /iK < ATksz < 30 nK, and -148 /xK < AT < 110 /iK, 
respectively. 



Figure|2]shows the convergence profiles k{9) averaged over 
annulus as a function of angular separation from the cluster 
center Each line represents k{9) from the average conver- 
gence fields obtained by stacking as many clusters indicated 
in the legend. Projected NFW profiles with truncation at i?vii 
provide a good approximation to the average cluster mass pro- 
file (thin solid) obtained by counting only the cluster members 
identified by the FoF algorithm. However, weak lensing mea- 
sures the projected mass distribution including contributions 
from interlopers that happen to lie between the lensing clus- 
ter and the observer The average convergence profile from 
the stacked field, therefore, provides the cluster-mass cross- 
correlation (thick solid). 



B. Secondary anisotropies by clusters 

On cluster scales, the CMB temperature and polarization 
fields can be approximated as large-scale gradient fields, and 
their gradient directions are nearly uncorrelated. Gravitational 
lensing by massive clusters induce dipolelike wiggles around 
the clusters to the CMB temperature and polarization gradi- 



ent fields, and these features are used to measure the lensing 
effect | |22[] . However, there exist several sources of contami- 
nation that mimic the lensing signature and complicate its re- 
construction process. 

Hot electrons in massive clusters scatter off CMB pho- 
tons, giving rise to a spectral distortion of the otherwise 
Planck distribution of the CMB temperature field. This ther- 
mal Sunyaev-Zel'dovich (tSZ) effect 1231 can be as large as 
AT ~ 500/iK and it is a major source of contamination, 
significantly exceeding the lensing signals from the clusters. 
However, with its distinctive spectral signature the tSZ ef- 
fect can be removed by multi-frequency observations and here 
we ignore the possible contamination of the residual tSZ ef- 
fect arising from an imperfect cleaning process. However, 
the scattering of hot electrons in massive clusters also gives 
rise to the Doppler effect due to the bulk motion of the clus- 
ters. Though smaller than the tSZ effect, this kinetic Sunyaev- 
Zel'dovich (kSZ) effect |24] is comparable to the lensing sig- 
nals ATksz — 30/iK and its identical spectral dependence 
makes it hard to separate from the intrinsic CMB temperature 
anisotropies or the lensing signature by clusters. 

In contrast to CMB temperature anisotropies, CMB polar- 
ization anisotropies have relatively little contamination from 
the tSZ and kSZ effects; There are a couple of contamination 
sources in CMB polarization anisotropies such as the scatter- 
ing of the kinetic quadrupole in the electron's rest frame, the 
scattering of the anisotropic CMB photons from the tSZ ef- 
fect, and double scattering within the cluster (e.g., ll25ll26ll ). 
However, compared to the lensing effect ^ 1/iK in polariza- 
tion, the effects of these contaminants are negligible <C 0. 1/iK 
as they are proportional to /3^t, yr, and /3t^ with the tangen- 
tial velocity /? — vt/c ^ 10"^ and y = rkTe/nieC^ ^ lO^'^ 
(e.g., I27ii28il ). Therefore, we only consider the kSZ effect as 
a contaminant in CMB temperature anisotropies and no con- 
taminant is assumed in CMB polarization anisotropies in the 
remainder of the paper. 

Figure [3] illustrates the kSZ effect around the massive clus- 
ter at z = 1 in Fig.[T] The top panels show the Thomson scat- 
tering optical depth r and the line-of-sight velocity field v\os 
of free electrons. We assume that the helium mass fraction is 
^Hc — 0.24 and that gas is fully ionized throughout the simu- 
lation box. While the gas distribution in massive clusters may 
have a velocity dispersion lower than the dark matter distri- 
bution, here we simply assume that the gas distribution traces 
the dark matter distribution in phase space with the universal 
mass fraction flb/ildm and investigate its impact in Sec. lIVBl 
Therefore, the top left panel appears identical to the conver- 
gence field (top left) in Fig. [l] as they are both proportional 
to the projected mass density. The cluster is moving toward 
the observer at the mean wios ~ —300 km s~^ and its uios dis- 
tribution is random with the rms a^ = 640 kms^ , though 
the average v\os seen in Fig. [3] is smaller due to the inclusion 
of random interlopers along the line-of-sight within the finite 
angular size of detector pixel. 

The bottom panels show the kSZ effect ATksz/T = 
-~T vios/c with and without the intrinsic CMB temperature 
field. Most of the intergalactic medium is transparent (r ~ 0) 
and there is little kSZ effect therein, except a few small blobs 
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FIG. 4: Lensing reconstruction of the convergence profile k{6) using improved quadratic estimators based on CMB temperature and polariza- 
tion fields in an ideal CMB experiment (apix — ^fwhm = 0). The initial model for the convergence field is assumed to have a uniform matter 
distribution k™(&) = 0.02 (dashed) and improved quadratic estimators are applied to 500 clusters, yielding an estimate of the convergence 
profile (thin solid) with its uncertainty (shaded) at each iteration, and refining the initial input model. All the improved quadratic estimators 
converge quickly to the true convergence profile (thick solid) in a few iterations. 



along the large-scale filaments. The kSZ effect ATksz — 
30/iK is highly concentrated at the cluster center and dimin- 
ishes fast as the matter density of clusters falls off rapidly. 
However, the contamination from the kSZ effect is apparent 
in the bottom right panel and its impact increases with redshift 
as clusters at high redshift are more compact and the line-of- 
sight velocity decreases only with (1 + z)^^/^. 



IV. RECONSTRUCTING CLUSTER-MASS 
CROSS-CORRELATION 



Using the numerical simulation in Sec. Hill we demonstrate 
the applicability of our improved quadratic estimators to real- 
istic clusters. In Sec. lIV Al we first test five improved quadratic 
estimators in an ideal CMB experiment, and we then compare 
their performance in realistic CMB experiments. In Sec. lIV Bl 
we discuss the impact of the kSZ contamination on the lensing 
reconstruction by using CMB temperature anisotropics and 
comment on a way to improve the reconstruction in the pres- 
ence of the kSZ effect. 



Improved quadratic estimators using CMB polarization 
anisotropics 



Here we test our improved quadratic estimators against the 
numerical simulation described in Sec. [Ill] Figure |4] shows 
the convergence profiles reconstructed by applying improved 
quadratic estimators (labeled in the legend of each panel) to 
500 clusters at z — 1 in an ideal experiment with no detector 
noise (iTpix = 0) and telescope beam (^fwhm = 0). Note 
that each cluster we stack has different shapes. We first ap- 
ply improved quadratic estimators by adopting an initial con- 
vergence model K"^{d) — 0.02 (horizontal dashed), i.e., no 
lensing signature in the observed CMB temperature and polar- 
ization anisotropics. Each estimator takes the assumed initial 
model and computes the change An from the assumed model, 
yielding an estimate (thin solid) of the convergence profile. 
The first estimates of each of the estimators from the uniform 
matter distribution are already close to the true convergence 
profile (thick solid). Having adopted the initial model with- 
out lensing effects, our improved quadratic estimators oper- 
ate as standard quadratic estimators, and they are biased low 
when the lensing effect is large near the massive clusters Ja- 
10], while the nonlinear effect is relatively mitigated for iTE, 
iTB, and iEB estimators, since the lensing contribution is par- 
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FIG. 5: Relative performance comparison of the improved quadratic 
estimators for reconstructing the convergence profile. The left panel 
shows the statistical uncertainties in the mean convergence profile 
obtained by applying each estimator to 1000 clusters in the fiducial 
CMB experiment with (jpix = 3/iK and Spwhm = 0.'5. The right 
panel shows the degradation in the total signal-to-noise ratio of each 
estimator as a function of the detector noise, compared to that of the 
iTT estimator in the fiducial experiment. 



tially canceled due to the oscillating nature of the cross power 
spectrum Cf and there is no B-mode polarization anisotropy 

cp = 0. 

Next we take the estimate (thin line) of the convergence 
profile in the first trial as our initial convergence model for 
the next iteration, and we repeat the iteration process until the 
estimates converge. In just a few iterations, the initial con- 
vergence model of the uniform matter distribution is quickly 
reshaped and all the estimators in each panel converge to 
the true convergence profile without any detectable bias be- 
yond the cluster virial radius (vertical dotted). The uncertain- 
ties (shaded region) in the mean estimate become smaller as 
the assumed convergence model is refined at each iteration; 
their minimum is achieved when the estimates converge to the 
true convergence profile and it is set by the intrinsic fluctua- 
tions of CMB temperature and polarization anisotropies. For 
a faster convergence, one can start with a more realistic ini- 
tial model of the convergence field motivated by other obser- 
vations rather than the uniform matter distribution assumed 
here. Note that while we characterize our reconstruction for 
the cluster-mass cross-correlation in terms of the convergence 
profile k{9) averaged over the annulus, it is the 2D conver- 
gence field K(n) that we reconstruct using improved quadratic 
estimators (similarly for the initial models). 

Figure |5] compares the relative performance of each im- 
proved quadratic estimator in reconstructing the convergence 
profile in the fiducial CMB experiment with Upix — 3/^K and 
^FWHM = 0.'5 and quantifies its signal-to-noise ratio for dif- 
ferent CMB experiments with varying levels of detector noise. 
The left panel plots the statistical uncertainty in the mean esti- 
mate of the convergence profile for each estimator from 1000 
clusters, after a few iterations when the estimates have con- 
verged. The statistical uncertainties are averaged over the log- 
arithmic radial bin, comparing the relative performance and 

1/2 

they decrease in proportion to N^{ with A^ci being the num- 
ber of clusters stacked. 

The performance difference of each estimator in Fig. |5^ 
originates from their noise power spectrum in Eq. (|9]l and the 



noise power spectrum TV; ' of each estimator, or equiv- 
alently the normalization N[^^ in Eq. (fT2l l. is determined 
by the intrinsic CMB power spectra C/^^, C{^ , and C^, 
given experimental specifications (detector noise and tele- 
scope beam size). Especially, the relevant information about 
the lensing reconstruction is contained around I ~ 1/6'ci ^ 
3000. Since the intrinsic CMB power spectra decay expo- 
nentially, the noise power spectrum N^ ' for iTT estima- 
tors at I ^ l/^ci is set by the sum of the ratios Cf /C^ ' 
at / < Icrit and /crit is the scale at which the ratio becomes 
unity (1/6'ci < ^ctr)- Analogously, the relevant ratios can be 
obtained by using Eq. ( fT2b as C/^ /C^ ' for iEB and iEE es- 
timators and {Cf)'^/Ci ' C^ ' for iTB and TE estimators, 
though the critical scale Zciit is similar to the cluster scale, 
1/6'ci — ^ciit for the polarization estimators. Since there is 
no intrinsic B-mode polarization anisotropy, the noise power 
spectra of iEB and iTB estimators are somewhat smaller than 
those of iEE and ITE estimators on cluster scales. Therefore, 
the lensing reconstruction of the cluster convergence profile 
can be best achieved by iTT estimators, and iEB, iEE, iTB, 
and iTE estimators have larger variance in the sequential or- 
der for the fiducial CMB experiment as shown in Fig.|5^. 

To quantify the signal-to-noise ratio of the improved 
quadratic estimators and their dependence on the detector 
noise, we compute the covariance matrix of each estimator 
as 

Cxy{0, 0') = ( [^""^{0) - ^{0)] [^^^(0') - n{e')] ) , 

(14) 
and the signal-to-noise ratio for each experiment specifica- 
tions is then 

^xxY = E ^(^) [^^^(^' ^')r' <(^') ■ (15) 

e,e' 

Since the lensing effect on CMB power spectra is to convolve 
them with its potential power spectrum C[, the lensing esti- 
mators are intrinsically non-local and the covariance matrix 
is non-diagonal 1 10]. The covariance matrix is computed by 
averaging over 50,000 clusters to guarantee its numerical con- 
vergence. 

Figure [Sj) plots the ratio of Ax 
the fiducial experiment to Ax^y of each estimator as a func- 
tion of detector noise. As in Fig. |5^, the general trend of 
the performance of each estimator remains unchanged over 
a wide range of detector noise; as the detector noise increases, 
the performance of all the estimators is progressively de- 
graded. However, the signal-to-noise ratios dramatically im- 
prove at (Tpix < 3/iK for the estimators using CMB polariza- 
tion anisotropies, since at the cluster scale the intrinsic CMB 
polarization anisotropies are small and easily dominated by 
the detector noise. 



B. Kinetic Sunyaev-Zel'dovich effect 



'^rp of the TT estimator in 



Now we use the simulation in Sec. IIIIBI with its velocity 
distribution to investigate the effect of the kSZ contamination 
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angular separation (arcminute) 

FIG. 6: Impact of the kinetic Sunyaev-Zel'dovich (kSZ) effect on the 
lensing reconstruction. For the fiducial experiment with Sfwhm ~ 
O.'S and (jpix = 3 /^K, we use the simulation in Sec.|lIl]to generate 
10,000 massive clusters at z = 1 and stack the convergence profiles 
reconstructed by applying iTT and iEE estimators. Left panel: As- 
suming that the gas distribution traces the dark matter distribution 
in phase space from the simulation, we compute the kSZ contami- 
nation as a function of the ionized gas fraction F in Eq. l ll6t . The 
reconstructed convergence profile is progressively biased in the pres- 
ence of the kSZ contamination in CMB temperature anisotropics. 
Right panel: Convergence profiles with iEE estimators as a function 
of telescope beam. No contamination from the clusters is assumed in 
CMB polarization anisotropies. Reconstructed convergence profiles 
are shown as dashed lines with their uncertainty (shaded), while the 
true convergence profile is shown as solid lines. 



illustrated in Fig. [3] on the lensing reconstruction. Note that 
we consider only the kSZ effect as a possible contamination of 
CMB temperature anisotropies and no other contamination is 
assumed for CMB polarization anisotropies, except the lens- 
ing effect itself (see, Sec. lIIIBl for estimates of possible con- 
tamination in CMB polarization anisotropies). 

In massive clusters, baryons make up a universal mass frac- 
tion J7b/51.,„ — 0.138 and they exist predominantly in the form 
of ionized gas in the intracluster medium. However, stars and 
galaxies in massive clusters contain a non-negligible fraction 
of baryons, and X-ray observations show that the gas mass 
fraction is /gas = 0.117 within r25oo, which corresponds to 
~ 0.1i?vir: Approximately 85% of the baryons in massive 
clusters contributes to the kSZ effect (see, e.g., l29l l30ll for 
estimating the gas mass fraction). Here we simply adopt the 
ionized gas fraction F as a free parameter and the free electron 
number density is then obtained from the simulation by 



Fx 






dm 



mdn 



1- 



njc 



(16) 



where nip is the proton mass. Since the kSZ effect is propor- 
tional to the product of the free electron number density and 
the Une-of-sight velocity, a lower value of F can represent a 
lower velocity distribution of baryons in massive clusters with 
the universal mass fraction ilb/flm- 

Figure |6^ shows the impact of the kSZ contamination on 
the lensing reconstruction (dashed) from iTT estimators in 
the fiducial experiment. When most baryons are contained 
in the ionized gas (F = 1 and F = 0.85), the kSZ contam- 
ination is ATksz — 30/iK and 25/iK at the cluster center. 



respectively, and its effect is substantially larger than the typ- 
ical change AT ~ 10/^K arising from the lensing effect by 
clusters. Consequently, the lensing reconstruction is signifi- 
cantly biased in the presence of the kSZ contamination. For 
the other two cases, in which large fraction of the baryons are 
in stars (F — 0.5 and F = 0.25), the kSZ contamination 
ATksz ~ 15/iK and 7.5/iK is smaller and its impact is re- 
duced, though the bias in the lensing reconstruction prevails 
over a range of separations. 

Compared to CMB temperature anisotropies, CMB polar- 
ization anisotropies are relatively free from secondary con- 
tamination and the lensing effect is the dominant source of 
secondary anisotropies. However, the current CMB experi- 
ments such as SPT and ACT lack the ability to measure polar- 
ization anisotropies, and other upcoming CMB experiments 
may not be optimized for measuring the cluster lensing sig- 
nature in polarization anisotropies. Figure |6}5 illustrates the 
lensing reconstruction by using iEE estimators in future CMB 
polarization experiments with larger telescope beam size (B- 
mode polarization measurements are even harder due to the 
vanishing signal). No significant bias develops for iEE esti- 
mators within the virial radius. However, as the angular res- 
olution decreases, small scale structure is smoothed and the 
reconstructed convergence profile is the true convergence pro- 
file convolved with the telescope beam (dashed line). 

Finally, with the current CMB experiments capable of mea- 
suring only temperature anisotropies on cluster scales, we 
consider a way to mitigate the impact of the kSZ contami- 
nation in the lensing reconstruction using CMB temperature 
anisotropies. Since the kSZ effect arises from the Thomson 
scattering off free electrons in clusters, its contamination is 
centrally concentrated as seen in Fig. [3] scaling in proportion 
to k{9). The lensing effect, however, extends well beyond the 
innermost region of clusters; for example, the kSZ effect falls 
off as 1/9 in singular isothermal clusters, while the lensing 
effect remains constant throughout the cluster region. There- 
fore, the simplest way to reduce the kSZ contamination is to 
suppress or purge the central region of clusters, where the kSZ 
effect is strongest (see, e.g., JSlll ). Since the lensing estima- 
tors are non-local and the information on the central region 
is shared by lensed CMB temperature anisotropies over the 
cluster virial radius lUOtl . the lensing reconstruction is possi- 
ble even with the pixels at the innermost region of clusters 
masked out. However, with our quadratic estimators built in 
Fourier space, translational invariance prevents from discrim- 
inating the central region. We leave the further investigation 
of how masking the central region affects the performance of 
our estimator for future work. 



V. DISCUSSION 

We have generalized the maximum likelihood estimators 
flOll for reconstructing the cluster-mass cross-correlation to 
utilize delensed CMB polarization fields, and we have tested 
our new maximum likelihood estimators for the lensing recon- 
struction against numerical simulations. The observed CMB 
temperature and polarization fields are delensed based on an 
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initial model of the convergence field and our maximum like- 
lihood estimators (called improved quadratic estimators) pro- 
vide improvements over the assumed initial model. This pro- 
cess can be iterated until numerical convergence is achieved. 
Compared to the standard and modified quadratic estimators 
DMIMI, our improved quadratic estimators have no free param- 
eters and provide unbiased reconstructions even in the regime 
where the linear approximation in the lensing effect breaks 
down. 

We have adopted smoothed hydrodynamic simulations to 
model realistic clusters, and the improved quadratic estima- 
tors that use six different combinations of CMB temperature 
and polarization fields can reconstruct the underlying mat- 
ter distribution in a non-parametric way to obtain the cluster- 
mass cross-correlation. Their ability to reconstruct the cluster- 
mass cross-correlation is determined by the ratio of the intrin- 
sic CMB power spectrum to the detector noise power spec- 
trum on cluster scales. Given the same experimental spec- 
ifications for measuring CMB temperature and polarization 
anisotropics, iTT estimators can reconstruct the cluster-mass 
cross-correlation by a factor of 3 in the signal-to-noise ratio, 
better than iEB estimators, which are better suited for recon- 
structing large-scale structure. 

For the gas mass fraction /gas — 0.117, the numerical 
simulations show that the kinetic Sunyaev-Zel'dovich (kSZ) 
contamination is ATksz — 30/iK at the innermost region of 
clusters. As clusters are more compact at high redshift with 
weak redshift dependence of their peculiar velocity, the spec- 



trally indistinguishable kSZ contamination in CMB tempera- 
ture anisotropics poses a significant challenge to the lensing 
reconstruction using iTT estimators. However, polarization 
estimators such as iEB estimators or more practically iEE es- 
timators can be used for the cross-check of the lensing recon- 
struction. Especially, these polarization estimators are more 
desirable given that there is relatively little kSZ contamina- 
tion in CMB polarization anisotropics and they perform better 
than iTT estimators in an ideal CMB experiment with no de- 
tector noise and telescope beam. However, with the prospect 
of CMB polarization experiments competitive with the current 
CMB temperature experiments far in the future, masking the 
central region of clusters to mitigate the kSZ contamination 
will be needed to take full advantage of precise CMB temper- 
ature measurements in the lensing reconstruction. 
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